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By  referencing  the  Lettan  Method  ,  the  verti¬ 
cal  distribution  of  turbulent  parameters  -  momentum 
transfer  coefficient  (k^,  friction  velocity  u^*, 
turbulent  shear  stress  ^  and  mixing  length  1,  were 
calculated  from  the  observed  data  of  average  wind 
profiles  in  the  neutral  planetary  boundary  layer 
/PBLt.  The  results  were  compared  with  research 
results  from  abroad  and  they  were  basically  in  fair 
agreement.  Some  of  the  parameters  should  be  modi¬ 
fied. 

/r 

I .  PREFACE 

The  basic  characteristic  of  the  planetary  boundary  layer 
is  its  turbulency  in  movement.  Turbulent  motion  is  resulted 
from  the  combined  effects  of  "vortex"  groups,  which  span  quite 
wide  in  the  time  and  space  scales,  under  complicated  boundary 
conditions.  To  directly  measure  the  turbulent  flow  field 
requires  high-tech  devices  and  is  generally  very  difficult. 
Therefore,  the  development  and  application  of  turbulent  sta¬ 
tistical  theories  has  been  limited.  Moreover,  using  the  K 
theory  to  solve  the  closed  equations  is  also  not  very  easy. 
Another  method  is  to  use  the  observed  average  wind  and  temper¬ 
ature  profiles  to  calculate  the  distribution  of  turbulent 
characteristics  under  assumed  conditions.  This  is  valuable 
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for  both  theoretical  and  experimental  studies. 


This  paper  used  the  field  wind  data  measured  by  the  bivane 
anemometer  at  Yizheng  of  Jiangsu  Province  on  April  29  and 
August  12,  1982  to  calculate  the  vertical  distribution  of 
turbulent  transfer  coefficient  ,  mixing  length  1,  friction 
velocity  u*  and  turbulent  shear  stress  T  in  the  neutral 
planetary  full  boundary  layer.  It  provided  a  base  for  the 
research  of  turbulent  structures  in  the  planetary  boundary  layer. 


II.  CALCULATION  METHOD 

We  referenced  the  Lettau  Method  Assuming  the  atmosperic- 

layer  structure  is  neutral  and  with  positive  pressure,  and  the 
horizontal  flow  field  is  steady  with  uniform  linear  motion, 
then  the  motion  equations  can  be  simplified  to: 


A. 


—  p/(t  —  vt) 
#»/(•*  ~  "t ) 


(1) 


where  and  ty  are  the  horizontal  components  of  turbulent  shear 
stress  respectively,  u  and  V  are  the  horizontal  components 
of  the  average  wind  velocity  field,  u  and  V  the  horizontal 
components  of  ground  rotation  wind.  f  is  the  Clarke  parameter 
and  f3  is  the  density  of  air.  Based  on  semi-experience  theory: 


T  =  Am 


av 

dz 


that  is, 


(2) 


Substituting  Equation  (2)  into  Equation  (1)  and  not  considering 
the  vertical  changes  of  density,  then  the  integration  with  respect 
to  alonq  the  ground  rotation  wind  direction  gives: 
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Take  two  special  layers  at  altitudes  2,  and  ^  >  which  are 
the  altitudes  where  u  and  V  reach  their  maximum  values,  then 
the  components  of  ground  shear  stress  can  be  written  as: 


,  i*1' 

r„  **■  PI  1  t/ds 


0> 


How  is  Ug  determined?  If  the  angle  0{Q  between  ground  rotational 
wind  and  ground  wind  is  given,  and  assuming  the  direction  of 
shear  stress  very  close  the  ground  is  the  same  as  that  of  the 
ground  wind,  then 

fjo  =  r«tancb 

Substituting  Equation  (5)  into  the  second  equation  of  (4) , 
we  get 

u,  “  —  udz  +  una,  |  vdz\  '  (6) 


Now  the  main  problem  is  the  determination  of  <xe  .  Trial  and 
error  method  can  be  used.  Assuming  turbulent  transfer  coeffi 
cient  is  a  standard  quantity,  i.e. 


*,•*>=*  *« 


(*) 


where  1  is  the  mixing  length.  An  initial  c*0  value  is  first 

given.  The  observed  wind  velocity  can  be  divided  into  u(z  ) 

and  "V (  z ) ,  and  z,  and  z~  are  located,  u  is  calculated  from 

1  z  g 

Equation  (6) .  T» o  can  be  obtained  by  substituting  ug  into 
Equation  (4),  and  TV©  obtained  from  Equation  (5).  Tx  (z) 
and  7*y  (z)  can  then  be  obtained  from  Equation  (3)  ,  and  finally 
kx  and  k^  are  calculated  from  Equation  (2).^ovalue  is  con- 
tinously  changed  until  (kx-k^)M  — *  0  (the  subscript  M  repre¬ 
sents  average  value  with  respect  to  height) .  Thus  accurate 
°<o  is  obtained.  Meanwhile,  the  turbulent  shear  stress 
T  (Tx,  Ty)  at  various  altitudes  is  obtained,  and  so  is  the 
turbulent  transfer  coefficient  Jcm .  Then  mixing  length  can  be 
solved  from 


-» -  t«VK£/ + (£  yrr 


(-■) 


The  friction  velocity  can  be  calculated  according  to  its  defi¬ 
nition:  !<•(*)  *=•  |rW/p]“  (»; 


where 


r(r)  -  [rlCO  +  rJCOl* 


For  details  of  the  calculations ,  refer  to  Reference  [1]. 

We  employed^&f gital  calculation  method,  and  the  Simpson 
formula  is  used  for  integration 

•i 

where  ,  b-a 

h= - 

m 

In  this  paper,  h=10m. 


Differential  calculation:  The  second-order  polynomial 
approximate  fitting  method  is  used.  The  coefficients  of  said 
polynomial  are  determined  by  the  least  mean  square  method 
using  the  observed  data  of  five  points.  At  a  straighter 
portion  of  the  profile,  simple  average  difference  formula  is 
used,  i.e.  <fy,  _  y-,.,  -  >,„  (J0) 

dx,  2  h 

We  improved  the  method  for  determining  ground  wind  direction. 
Instead  of  extrapolating  from  wind  data  measured  by  the  small 
balloon,  accurate  anemometer  for  direct  measurement  of  instan¬ 
taneous  wind  velocity  and  wind  direction  was  used.  The  above 
calculations  were  conducted  in  a  PC-1500  computer. 


III.  OBSERVATION  BACKGROUND  AND  MACROSCOPIC  METEOROLOGICAL 
CONDITIONS 


The  paper  primarily  used  the  wind  velocity  profile  data 
measured  by  the  bivane  anemometer  at  Yizheng  of  Jiangsu  Pro¬ 
vince  on  August  12,  1982.  The  terrain  surrounding  the  in¬ 
fluence  by  the  Chang  Jiang,  data  leaning  toward  north  wind 
were  selected. 

From  8:30-17:30  on  August  12,  there  were  12  groups  of 
wind  velocity  profiles.  It  could  be  seen  from  the  meteorological 


diagram  that  the  isobars  on  that  day  in  said  area  were  quite 
straight  and  that  the  isotherms  were  basically  parallel.  There 
were  no  changes  in  the  meteorological  system,  and  the  flow  field 
was  rather  stable.  At  14:00  the  low  altitude  air  observations 
gave:  average  vertical  temperature  decreasing  rate  was 

-0.61°C/100m  under  1  kilometer,  average  temperature  gradient 
was  0.43°C/100m.  Thus  the  atmospheric  structure  was  neutral. 

The  observation  field  and  macroscopic  meteorology  matched  the 
aforementioned  assumed  conditions. 


IV.  ANALYSIS  OF  CALCULATED  RESULTS 


First  of  all,  the  trial  and  error  method  was  used  to  de¬ 
termine  Cta  value.  On  August  12,  ground  barometric  pressure 
P„=976.2  milibar,  temperature  Tn=26.5°C;  calculated  density 

p  =1.24  x  10  g/cm  ;  latitude  0  =32°;  the  Clarke  paramter 
-4  -1 

f=0.773  x  10  sec  .  Table  1  gives  calculated  y results  for 
different  Oc  values. 


Table  1 
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Key:  (a)  degree;  (b)  meter;  (c)  dynes/cm  . 

(kx~ky)M  in  the  table  represents  the  average  transfer  co¬ 
efficient  under  the  altitude  of  400  meters.  It  can  be  seen 
from  the  above  table  that  the  angle  between  ground  rota¬ 
tional  wind  and  ground  wind  should  be  selected  as  24.4°  with  a 
ground  rotational  wind  velocity  of  14.71m/sec. 
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The  August  12  spirals  diagram  of  the  average  wind  velocity 

profiles  is  as  shown  in  Figure  1.  The  profile  of  the  wind 

velocity  component  is  as  shown  in  Figure  2.  In  Figure  2,  the 

Y-axis  is  the  net  height.  Compared  with  Lepzig  data,  the 

trend  of  the  two  distributions  is  similar.  Figure  3  shows  the 

spirals  of  turbulent  stress  calculated  from  the  average  wind 

velocity  profile  and  (A.,  and  they  are  in  fair  agreement  with 

[11  u 

those  of  Lettau 

Using  the  aforementioned  method,  the  vertical  distribution 

of  friction  velocity  u*,  turbulent  transfer  coefficient  k_ 

1 

and  mixing  length  were  calculated. 


Fig.  1.  Spiral  of  the  observed  wind 
velocity  on  August  12,  1982,  with  V 
direction  as  X-axis  direction.  3 
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Fig.  2.  Observed  wind  velocity  component 
comparing  with  Lepzig  data 
Dotted  line:  Lepzig 

Solid  line:  observed  data  of  this  paper. 


1  dyne /cm 


Fig.  3.  Calculated  turbulent  shear 
stress  spiral  (<3Cq=24.4°) 

The  values  at  various  altitudes  from  the  ground  to  1 
kilometer  are  listed  in  Table  2  and  are  as  shown  in  Fig. 4. 

The  figure  gives  the  km  profiles  of  August  12  and  April  29 
and  those  calculated  from  the  Lepzig  data.  The  comparison 
shows  that  the  trend  of  their  distribution  is  in  fair  agree¬ 
ment  (especially  obvious  above  400  meters).  On  the  average, 
the  km  values  of  this  paper  are  greater  than  those  of  Lepzig' s. 
The  difference  mainly  occurs  at  a  lower  altitude.  The  Lepzig 
1%,  curve  increases  very  quickly  under  100  meters  with  the 
maximum  occurring  at  200  meters.  The  curve  of  our  measuring 
point  increases  slowly  under  100  meters,  then  increases  sharply 
reaching  its  maximum  at  around  250  meters.  The  curves  for 
April  29  and  August  12  match  well  under  200  meters  (the  dis¬ 
crepancy  at  a  high  altitude  was  due  to  the  flow  field  of 
April  29  not  meeting  the  assumed  conditions  fully,  and  effects 
of  thermal  wind  were  present  at  a  higher  altitude) .  Moreover, 
the  Lepzig  profiles  were  the  average  profiles  of  28  observed 
data  in  one  day,  yet  the  profiles  of  this  paper  for  August  12 
were  the  average  of  12  observed  data  on  that  day,  making  the 
curve  not  as  smooth.  The  general  trend,  however,  is  in  fair 
agreement.  Therefore,  the  Lettau  method  is  creditable.  It 
can  be  seen  from  comparing  Fig. 2  and  Fig.  4  that  the  distribu¬ 
tions  of  km  and  v  are  basically  similar.  Just  as  explained  by 
[11 

Lettau  ,  the  energy  of  turbulent  momentum  transfer  is  obtained 
from  the  "potential  energy"  of  horizontal  barometric  pressure 
through  the  effects  of  V*  . 


Table  2.  Vertical  distribution  of  wind  velocity 
component  and  km 
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The  vertical  distribution  of  the  calculated  mixing  length 
1  values  is  as  shown  in  Fig.  5.  From  the  ground  up,  1  increases 
with  altitude.  At  the  altitude  between  260-300m,  it  reaches 
the  maximum,  then  starts  to  decrease  slightly.  This  is  con- 
sistent-with  the  observation  of  RossbyL  .  Figure  5  also  gives 
the  calculated  values  based  on  the  Blackader^  and  Lettau  ^ 
experience  models. 
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Fig.  4.  Vertical  distribution  of 

km.—  .  August  12 - April  29  —•  — 

Lepzig. 


Fig.  5.  Vertical  distribu 
tion  of  1  .August  12 
Lettau  - Blackader. 


Blackader  (1962) 


1  =  2.7  X  lQ~'p,/f 


Lettau  (1962) 


t  = - — - 

1  +  ( */Ov< 


where  Zm  is  the  altitude  where  1  reaches  its  maximum,  and  x 
is  the  Carmon  constant. 


It  can  be  seen  from  Fig.  5  that  there  is  a  difference 
between  our  calculated  1  values  and  those  of  the  above  two  for 
mulae.  The  main  discrepancy  occurs  in  the  upper  portion  of 
the  boundary  layer  where  the  trend  of  1  distribution  is  quite 
different.  Formula  (11)  gives  1  values  which  increase  quickly 
from  the  ground  up  to  about  300  m,  then  the  increase  slows 
down  as  altitude  increases,  formula  (12)  gives  1  values  which 
increases  sharply  up  to  400m,  and  from  then  on  the  values  stay 
essentailly  unchanged.  The  1  values  of  this  paper  which  were 
calculated  from  the  observed  data  increase  from  the  ground  up. 
and  the  maximum  value  occurs  at  around  300m.  From  then  on, 
they  slowly  decrease.  If  formula  (12)  is  modified  as 

/  = - if -  (23) 

I  4  1 

and  its  x  is  set  as  0.35  which  was  proposed  by  Bussinger 
(1971),  then  the  calculated  values  of  formula  (13)  are  in 
better  agreement  with  those  of  the  observed  data  as  shown  by 
the  solid  line  in  Fig.  5. 


The  calculated  results  of  friction  velocity  u*  are  as 
shown  in  Fig.  6.  From  the  ground  up,  u*  decreases  linearly 
with  respect  to  altitude.  .The  experience  formula  for  u*  in 
the  neutral  tower  layer  proposed  by  Clarke ^ (1970)  is: 

■<*  —  —  6f:  (14) 

Based  on  our  calculated  results,  at  altitude  under  200m,  the 
following  formula 


-  5 lz 


will  fit  better.  As  for  the  basic  trend  of  linear  decrease  with 
respect  to  altitude  for  u*,  both  formulae  are  in  total  agreement 
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Fig.  6.  Vertical  distribution  of  friction 
velocity  u*  -  results  of  this  paper  - 
Clarke. 

V.  CONCLUSIONS 

(1)  Our  calculated  results  have  shown:  by  using  the 
Lettau  method,  it  is  feasible  to  calculate  turbulent  parameters 
of  the  neutral,  horizontally  uniform  and  steady  flow  field  from 
average  wind  velocity  profiles.  More  satisfactory  results  can 
be  obtained  as  long  as  observed  data  meeting  the  assumed  con¬ 
ditions  are  available. 

(2)  The  primary  results  of  the  vertical  distribution  of 
the  neutral  boundary  layer  turbulent  parameters  calculated  from 
observed  wind  velocity  data  are  as  follows: 

The  km  values,  from  the  ground  up,  increase  linearly  with 
respect  to  altitude.  The  maximum  values  occur  around  250m, 


decrease  gradually. 


The  1  values,  from  the  ground  up,  increase  linearly  to  the 
range  of  250-300m  where  the  maximum  occurs,  then  decrease 
gradually. 

The  u*  values  decrease  linearly  from  the  ground  up.  The 
rate  of  direct  decrease  is  5f. 

(3)  Our  calculation  results  are  compared  with  the  research 
results  from  abroad,  and  the  distribution  profiles  of  7",  u*, 

km  and  1  in  the  entire  boundary  layer  are  in  fair  agreement. 

But  there  are  discrepancies  in  specific  distribution  profiles, 
which  are  shown  when  comparing  formulae  (11)  and  (12)  with 
formula  (13);  formula  (14)  with  formula  (15);  and  km  profile 
in  Fig.  -4.  As  stated  before,  the  discrepancies  are  caused  by 
difference  in  ground  roughness.  In  addition,  our  observed 
data  for  August  12  were  obtained  from  12  groups  of  profiles  for 
the  period  between  08:30  and  17:30.  There  was  a  little  change 
in  the  temperature  structure  during  this  period.  By  analyzing 
the  synchromized  data  from  low  altitude  air  observations: 
0.4-0.8(  C/lOOm) ;  =0 . 45°C/100m.  They  did  not  satisfy  the 

neutral  conditions  completely. 

(4)  There  are  given  limitations  in  the  use  of  wind  data 
measured  by  the  bivane  anemometer  ~to  calculate  turbulent  para¬ 
meters: 

(a)  It  can  not  accurately  calculate  the  average  wind 
velocity  at  different  altitudes  for  the  same  time  instant  from 
the  trace  of  the  small  balloon  which  rises  at  a  constant  speed. 
First  of  all,  the  timing  of  measurement  for  different  altitudes 
can  not  be  strictly  "synchronized".  For  a  small  balloon  with 
a  rising  speed  of  lOOm/min,  the  time  lag  can  reach  15  minutes 
between  the  first  measurement  and  that  at  an  altitude  of  1.5km. 
Secondly,  the  time  intervals  of  the  average  values  are  not  very 


representative.  If  more  measuring  points  are  added,  then  the 
time  interval  for  taking  readings  must  be  shortened (i.e. the 
average  time  interval) .  Our  observation  took  4  readings  in 
one  minute,  i.e. every  wind  velocity  value  represents  the 
average  wind  velocity  within  a  25m  altitude  range.  Obviously 
this  time  interval  is  quite  small  and  can  not  include  turbu¬ 
lent  vortexes  of  a  wider  scale.  Therefore,  the  profiles  are 
not  smooth  enough  and  must  conduct  numerous  repeat  observations 
to  obtain  the  average.  Yet  it  is  difficult  to  maintain  stable 
temperature  structure  and  macroscopic  meteorological  conditions 
during  the.  repear  observation  period,  thus  causing  the  applica¬ 
tion  of  said  method  to  have  given  limitations. 

(b)  The  baseline  selected  for  bivane  anemometer  observa¬ 
tions  must  be  perpendicular  to  wind  direction.  Wind  direction 
close  to-  the  ground  and  the  boundary  layer  usually  would  shift 
20°-30°  or  more  as  altitude  increases.  Therefore,  the  same 
baseline,  which  is  proper  for  use  on  the  ground,  becomes  not 
perpendicular  to  wind  direction  as  the  balloon  rises  to  the 
top  portion  of  the  boundary  layer  and  the  wind  direction  shifts. 
Even  a  difference  of  just  a  few  degrees  in  the  two  bivane 
anemometer  direction  readings  will  cause  great  calculation 
errors,  and  the  vertical  projection  method  or  vector  method 
must  be  used.  If  the  directional  angles  differ  little  and 
so  do  the  attack  angles,  then  the  above  two  methods  will  do 
little  good.  We  imagine,  if  two  baselines  are  used  simultan¬ 
eously  along  with  a  trivane  anemometer  to  conduct  observations, 
the  aforementioned  defects  can  expect  to  be  alleviated,  thereby 
raising  the  accuracy  of  wind  velocity  profile  measurements  in 
the  boundary  layer. 

The  observed  wind  velocity  profile  data  used  in  this  paper 
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